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Abstract. We prove a theorem on algebraic osculation and we apply our 
result to the Computer Algebra problem of polynomial factorization. We con- 
sider X a smooth completion of and D an effective divisor with support 
dX = X\C? . Our main result gives explicit conditions equivalent to that a 
given Cartier divisor on the subscheme {\D\, Ob) extends to X. These oscula- 
tion criterions are expressed with residues. We derive from this result a toric 
Hensel lifting which permits to compute the absolute factorization of a bivari- 
ate polynomial by taking in account the geometry of its Newton polytope. In 
particular, we reduce the number of possible recombinations when compared 
to the Galligo-Rupprecht algorithm. 



1. Introduction 

This article is originally motivated by the wellknown Computer Algebra problem 
of polynomial factorization. We introduce here a new method to compute the 
absolute factorization of sparse bivariate polynomials, which is based on criterions 
for algebraic osculation in toric varieties. 

Since the eightees, several deterministic or probabilistic algorithms have been 
obtained to compute the irreducible absolute factorization of dense bivariate poly- 
nomials defined over a number field if C C. We refer the reader to [T7] and [5] 
for a large overview of the subject. In many cases, these algorithms are based 
on a Lifting and Recombination scheme - refered here as LR-algorithms - which 
detects the irreducible absolute factors of a polynomial / G K[ti,t2\ in its formal 
decomposition in .ftr[[ii]][t2]- Although this approach a priori necessites an expo- 
nential number of possible recombinations, people succeded in the last decade to 
develop LR-algorithms running now in a quasi-optimal polynomial complexity (see 
for instance [8], [1], [9] and the reference within). 

In general, a LR-algorithm first performs a generic afhne change of coordinates. 
When / has many zero coefficients in its dense degree d monomial expansion - we 
say that / is sparse - this step loses crucial information. In this article, we propose 
a new method which avoids this generical choice of coordinates. Roughly speaking, 
we obtain a toric version of the Hensel lifting process (in the spirit of [1]) which 
detects and computes the irreducible absolute factors of / by taking in account 
the geometry of its Newton polytope Nf. This permits in particular to reduce 
the number of possible recombinations when compared to the Galligo-Rupprecht 
algorithm [T5j . 

Let us present our main results. 

Algebraic Osculation. A natural way to take advantage of the Newton polytope 
information is to embed the complex curve of / in a smooth toric compactification 
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X of the complex plane. For a well chosen X, we can recover Nf (up to translation) 
from the Picard class of the Zariski closure C C X oi the afhne curve of / and we 
want to use this information. The boundary dX = X \<C^ of X is a. normal crossing 
divisor whose Picard group satisfies 

Pic{dX) ~ Pic(X). 

Thus, it's natural to pay attention to the restriction of C to some effective Cartier 
divisor D (more precisely, to the subscheme {\D\, Od)) with support \dX\. In order 
to detect the irreducible components of C, we need to find conditions for when a 
Cartier divisor on D extends to X. This is achieved in our main 

Theorem 1. Let X be a smooth projective compactification o/C^, whose boundary 
dX = X \ is a normal crossing divisor. Let D be an effective Cartier divisor 
with support \dX\ and let il^{D) be the sheaf of meromorphic 2-forms with polar 
locus bounded by D. 

There exists a pairing (•, •) between the group of Cartier divisors on D and the 
vector space H'^{X, il'j^(D)) with the property that a Cartier divisor j on D extends 
to a Cartier divisor E on X if and only if 

{^,^) = eH^{X,nj,{D)). 

The divisor E is unique up to rational equivalence. 

Not surprisingly, we'll construct such a pairing by using Grothendieck residues 
(see for instance ^1] where the authors study the interplay between residues and 
zero-dimensional subschemes extension). When X is a toric surface we obtain an 
explicit formula for (7, ^I^), generalizing a theorem of Wood [34 . To prove Theorem 
1, we compute the cohomological obstruction to extend line bundles from D to X 
and then we use the Dolbeault 9-resolution and residue currents to make explicit 
the conditions. 

Application to polynomial factorization. If two algebraic curves of fixed de- 
gree osculate each other with sufficiently big contact orders on some finite subset, 
they necessarily have a common component. This basic observation permits to 
derive from Theorem 1 an algorithm which computes the absolute factorizion of a 
bivariate polynomial /. The polynomial / is assumed to be defined over a subsfield 
if C C and we look for its irreducible decomposition over C. 

Under the assumption /(O) 7^ 0, we can associate to / a smooth toric completion 
X of whose boundary 

dX^Di + --- + Dr 

is a normal crossing toric divisor and such that the curve C <Z X oi f does not 
contain any torus fixed points of X. A Minkowski sum decomposition 

Nf=P + Q 

of the Newton polytope of / corresponds to a line bundle decomposition 

where Cp and Cq are both globally generated with at least one non trivial global 
section. 

We prove the following result (Subsection 3.2) 
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Theorem 2. There exists a unique effective divisor D with support \dX\ and 
rationally equivalent to C + dX . Let 7 be the restriction of C to D and suppose 
that Nf = P + Q. There exists an absolute factor q of f with Newton polytope Q if 
and only if there exists < 7' < 7 such that 

deg(7' • A) = deg/:Q|£,^, i = l,...,r 

and so that the osculation conditions hold for the pair (D,7'). The factor q is 
computed from 7' by solving a sparse linear system of 2Vol(Q) + cieg(7' ■ dX) 
equations and Card((5 H Z^) unknowns. 

When the facet polynomials of / are square free over K, we can derive from 
Theorem 2 a vanishing-sum LR-algorithm (Subsection 3.4). It first computes the 
Newton polytope decomposition 

Nf ^Qi + --- + Qs 

associated to the absolute decomposition of /. Then it computes the associated 
irreducible absolute factorization 

f = qix--- xqs 

with floatting calculous and with a given precision. The numerical part of our al- 
gorithm reduces to the absolute factorization of the univariate exterior facet poly- 
nomials and the algorithmic complexity depends now on the Newton polytope Nf 
instead of the degree d — deg(/). In particular, we fully profit from the combi- 
natoric restrictions imposed by Ostrowski's conditions Npq = Np + Ng (see |28j). 
This permits to reduce the number of possible recombinations when compared to 
the Galligo-Rupprecht algorithm p31. 

Finally, let us mention that Theorem 1 concerns also non toric C^-completions. 
In theory, this permits to exploit the information given by the non toric singularities 
of C along the boundary of X when / has non reduced facet polynomials (see 
Subsection 3.6). 

A formal study of algorithmic complexity, as well as questions of using non toric 
singularities information will be explored in a further work. 

Related results. Our method is inspired by an algorithm presented in [12] that 
uses generical toric interpolation criterions [32j in an open neighborhood of dX 
(see Subsection 3.5). By using combinatorics tools, the authors in [T| obtain a 
comparable Hensel lifting process which there too takes in account the geometry 
of the Newton polytope. Finally, let us mention [3], where the authors reduce the 
factorization of a sparse polynomial to smaller dense factorizations. 

Organization. The article is organized as follow. Section 2 is devoted to algebraic 
osculation. We introduce the problem in Subsection 2.1 and we construct the 
residue pairing in Subsection 2.2. We enounce precisely Theorem 1 in Subsection 
2.3 and we give the proof in Subsection 2.4. We give an explicit formula for the 
osculation criterions when X is a toric surface in Subsection 2.5. In Section 3, we 
pay attention to polynomial factorization. We prove Theorem 2 in Subsections 3.2 
and 3.3 and we develop the sketch of a sparse polynomial factorization algorithm 
in Subsection 3.4. We compare the underlying algortihmic complexity with related 
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results in Subsections 3.5 and discuss non toric information in Subsection 3.6. We 
conclude in the last Subsection 3.7. 

Aknowledgment. We would like to thanks Michel Brion, Stephane Druel, Jose 
Ignacio Burgos and Martin Sombra for their disponibility and helpfuU comments. 
We thanks Mohamed Elkadi and Andre Galligo who suggested us to pay attention 
to the interplay between toric geometry and sparse polynomial factorization. 

2. Algebraic osculation 

2.1. Notations and motivation. In all the sequel, {X,Ox) designs a smooth 
projective surface where 

A = Ao U \dX\ 

is the disjoint union of an affine surface Ao ~ with the support of a simple 
normal crossing divisor 

dX = Di + --- + Dr. 
We say that A is a completion of the complex plane with boundary dX. 

An osculation data on the boundary of A is a pair {D,j) where 

D = {ki + l)Di + --- + {kr + l)Dr 

is an effective divisor with support \D\ ~ \dX\ and 7 is a Cartier divisor on the 
subscheme (|-D|, Od)- By abuse of language, we will write D = (\D\,Od). 

An osculating divisor for (Z?, 7) is a Cartier divisor on A which restricts to 7 
on D. That is 

= 7, 

where i : — > A is the inclusion map. In other words, we are looking for a divisor E 
with prescribed restriction to the fc*'*-infinitesimal neighborhood of Di. In general, 
such an osculating divisor does not exist and we are interested here to determine 
the necessary extra conditions. 

We say that 7 has support [rj, where F designs the zero-cycle 7 • dX. Thus, 7 
can be uniquely written as a finite sum 

7 = ^ 7p, 

pe|r| 

each 7p being the restriction to D of a germ of an analytic cycle of A 

7p = div(/p) 

at p. When 7p is smooth and intersects transversally dX, a curve restricts to jp 
at p if and only if it has contact order at least kp with 7p, where fcj, + 1 is the 
multiplicity of D at p. This observation motivates the terminology of algebraic 
osculation. 

2.2. Residues. It's a classical fact that Grothendieck residues play a crucial role 
in osculation and interpolation problems. Let us mention for instance [12], 

[31] and [32] for interpolation results and [21] for the interplay between residues 
and subscheme extensions. Not surprisingly, residues will appear here too. 

Let il,x be the canonical bundle of A. We identify the line bundle il^{D) with 
the sheaf of meromorphic forms with polar locus bounded by D. Thus, any global 
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section 5* of rij^{D) restricts to a closed 2-form on Xq. Since Xq ~ is simply 
connected, there exists a rational 1-form on X such that 

For p G |r|, we denote by ifjp the germ of ip in the chosen local coordinates. Let 
hp be a local equation of D at p. Thus /ip-iAp is holomorphic at p. Suppose for a 
while that fp is holomorphic and irreducible. Then, following [50], we define the 
Grothendieck residue at p of the germ of meromorphic 2-form dfp A ipp/ fp as 



(1) resp 



dfp . 1 1 /" dfp 

— AV-p :=hni— — 2/ -^Aipp, 
fp n 1^0 (2i7r)2 /p 



where Mp(e) = {x close top, = ei, |/ip| = €2}. 

This definition does not depend on the choice of local coordinates [20 . By Stokes 
Theorem, it only depends on di/j — Moreover, the local duality Theorem [20] 
implies that (1) only depends on fp modulo (hp). That is, (1) depends on jp and 
not on the chosen lifting jp = {fp — 0}. By linearity, we can extend (1) to any 
germ of analytic cycle jp = div(/p) and it follows finally that (1) defines a bilinear 
operator 

' df 

(7, ^')p := reSp A^l^p 
L Jp 

between the group of Cartier divisor of D and the C- vector space H^{X, n\{D)). 
We refer to the proof of Theorem 1 (Subsection 2.4) for more details. 

2.3. Criterions for algebraic osculation. We keep previous notations. Our 
main result is the following 

Theorem 1. Let X he a smooth projective completion ofC'^ with a normal crossing 
boundary dX and consider an osculating data (D,^) on dX. 

1. There exists a Cartier divisor E on X which restricts to ^ on D if and only 

(2) (7' *)p = * ^ ^"(^' ^x{D))- 

p£\T\ 

The divisor E is unique up to rational equivalence. 

2. If moreover 7 is effective and 

H\X, Ox{E~D)) = Q, 
then there exists an effective divisor of X which restricts to 7 on D. 



The necessity of (2) follows from the Residue Theorem [20 . The difficult part 
consists to show that these conditions are also sufficient. The proof will be given 
in the next Subsection 2.4. Let us first illustrate Theorem 1 on a simple example. 

Example 1 (the Reiss relation). Let X = and consider a finite collection of d > 
smooth analytic germs 7p transversal to a line L C P^. Suppose that we look for 
an algebraic curve C C P^ of degree d which osculates each germ with a contact 
order > 2. This problem leads to the osculation data {D,j), where D = 3L and 7 
is the restriction to D of X]pe|r| 'yp- There is an isomorphism 

iJ ° (p2 , f^^^ (3L) ) ~ (p2 , Op2 ) ~ C 
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and we can check that a generator is given by the form ^! whose restriction to 
= \ L is equal to ^'|C2 = dti A dt2 in afhne coordinates (ti, t2)- Thus we can 
choose ip\c^ — tidt2. Letting ti = Ti/Tq and t2 = T2/T0 we obtain 

Ti{TodT2~T2dTo) 



n 



in the P'^ homogeneous coordinates [Tq : Ti : T2\- Up to an afhne change of 
coordinates, we can suppose that 7 is supported in the affine chart T2 7^ 0. In the 
new afhne coordinates x = T0/T2 and y = T1/T2, the hne L has local equation 
X = and = —ydx/x^. Moreover, we can choose a Weierstrass equation 

ip = y- (t)p{x) 

for the smooth germ 7p, where (pp G C{x}. By Cauchy formula, we obtain 



ydx A d(y - 4>p) 



x^{y - 4>p) 



reso 



dx 



1 



where reso is a univariate residue and <f>'p is the second derivative of <j)p. Finally, (2) 
is here equivalent to that 

(3) ^<^;'(o) = o. 

pe|r| 

For degree reasons, any osculating divisor is rationally equivalent to dL. Since 
H^{Op2[d — 3)) = 0, the relation (3) is finally equivalent to that there exists an 
osculating curve C for {D,j). 

It's easy to see directly the necessity of (3). An osculating curve is given by a 
degree d polynomial C — {f{x,y) — 0} that can be factorized 

f{x,y) = n " "pW) 
pe|r| 

in C{a;}[j/]. Since deg{f) = d, the sum X]pe|r| ''^pi^) ^ degree 1 polynomial and 
the relation 

pG|r| 

holds. If C osculates jp with contact oder 2, then (j)p and Up have the same Taylor 
expansion up to order 2, which directly shows necessity of (3). When we express the 
second derivative of Up in terms of the partial derivative of /, we recover the classical 
Reiss relation 18 . This result is obtained by Griffiths-Harris in [5D], Chapter 6. 

2.4. Proof of Theorem 1. The Cartier divisor 7 corresponds to a line bundle 
C G Pic(£') over D together with a global meromorphic section /. An osculating 
divisor for {D,j) corresponds to a line bundle C G Pic{X) together with a global 
meromorphic section / which restrict respectively to £ and / on D. 

The following lemma gives the cohomological obstruction for the extension of £. 
All sheaves are considered here as analytic sheaves and any sheaf on a subscheme 
Y G X is implicitly considered as a sheaf on X by zero extension. 

Lemma 1. There is a decomposition o/Pic(£') in a direct sum 

Fic{D) = Pic(X) e H\D, Od). 
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Proof. The classical exponential exact sequence exists for any curve of X (reduced 
or not, see p. 63). We obtain the commutative diagram 

^ Zx ^ Ox """"-^'^ o\ 

(4) ill 

Zd ^ Od 0*0^0 

where Zd C Od is the subsheaf of Z-valued functions (so that Zd — '^\d\) and 
vertical arrows are surjective restriction maps. Since X is rational, H'^(X, Ox) = 
for i > 0. Using the associated long exact cohomological sequences, we obtain the 
commutative diagram 

^ Pic(X) ^ H'^{X,Zx) ^ 

1 Ir ir' i 

HHD,Zd) ^ H\D,Od) ^ Vic{D) H H^{D,Zd) ^ H^{D,Od). 

Here r, r' are restriction maps, 5d^ 5x are the coboundary maps corresponding to 
Chern classes on D and X (see [5j , Ch. 1 for the non reduced case) and e is induced 
by the exponential map. 

We claim that r' is an isomorphism. Let j : X^ X he the inclusion map of 
Xo = X\\D\. The short exact sequence 

(5) 0^j,(Zxo) ^Zpi -.0 
gives rise to the long exact cohomological sequence 

Hl[Xo) H\X,Zx) ^ HWDIZ\d\) ^ Hl{X^) 

where i/*(Xo) is the compact support cohomology of Xq. The cohomology H*{Xq) 
is dual to the singular homology H^,{Xq) which vanishes in degree 1,2,3 since 
Xq ~ C'^, and the claim follows. Thus 

5b o r — r' o 5x 

is an isomorphism, and Pic(X) is a direct summand of Vic{D). The exponential 
cohomological sequence for X is exact in degree so that {X, Zx ) ^ {X, Ox ) 
is injective and H'^{X,Zx) — 0. We have just seen that H^{Xq) — 0, and finally 
(5) implies i?i(|D|,Z|c|) = 0. It follows that keriSo) =i H\D,Od)- □ 

Corollary 1. There is an isomorphism Pic(X) ~ Pic{dX). 

Proof. By the previous lemma, it's enough to show that H^{Dred, Oor-^a) — where 
dX = Dred is the reduced part of D. Since H^{Z,Zz) = for a finite set Z, there 
is a surjection 

r 

H\\D\,Z\d\)^^H\D,,Zd^). 

i=l 

Thus H^{Di,ZBi) — for all i and each Di is a rational curve. In particular, 
H^{Od,) = 0. Moreover, H^{Xo) = so that H°{X,Zx) ^ H°{\DIZ\d\) is 
surjective and \D\ is connected. This shows that |D| is a connected and simply 
connected tree of rational curves. There thus exists i so that Dred = Di + D', 
where n |Z3'| is a point and D' remains a connected and simply connected tree. 
By induction on the number of irreducible branches of Dred, it's enough to show 
that H^{Od,.^^) — H^{Od')- The short exact sequence 

^ Od„, ^ Od. ® Od' Oa^.jj' 
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gives the long exact sequence in coliomology 

^ H\Od^^,) ^ H\Od,)®H\D',Od') ^ H\Od,.d'). 

By assumption, Dred is a normal crossing divisor and Di ■ D' = {pt} as a reduced 
subscheme. The diagram then begins byO^C^C©C^C, which forces r to be 
injective. Since H^{OdJ = H^{0{pt}) = 0, this gives H\Od,,^) ^ H^{Od')- □ 

By the previous corollary, there exists a unique line bundle C £ Pic{X) so that 
C\QX — C.\dx s-nd L lifts to X if and only if 

Since ^^{Lq) = (by construction), there exists /3 e H^{D,Od) with e(/3) = £o- 
Moreover e is injective and 

Co^Od ^ [3 = 0. 

There remains to make explicit such a condition. We first use Cech cohomology. 
Suppose that Co is given by the 1-cocycle class g = {guv} S H^{D, CJ,), relative to 
some open covering U of X. We can suppose U fine enough to ensure a logarithmic 
determination 

huv ■.= :^logiguv) e OxiUnV), 

where guv G O^iU n V) is some local hfting of guv- Since S^ig) = 0, the classes 
huv & Od{U n V) of huv define a 1-cocycle class {huv} G H^{D,Od) which 
represents P (this definition does not depend on the choice of the liftings). Since 
X is rational, the structural sequence for D 

Ox{-D) -^Ox ^Od ^0 

gives rise to a coboundary isomorphism 5 : H^{D,Od) H^{X,Ox{—D)) while 
Serre Duality gives a non degenerate pairing 

H^{X, Oxi-D)) iJ°(X, nj,{D)) ^ H\X, ~ C, 

where Tr is the trace map [5,. We identify Ox{—D) with the sheaf of functions 
vanishing on D and (D) with the sheaf of meromorphic 2-forms with polar locus 
bounded by D. The Serre pairing ((5(/3), ^I^) with * e H"{X, il'^{D)) is represented 
by the 2-cocyle class 

C* ■= ^{huv + hvw + hwu)'^\unvnw^ G H'^{X, ^x) 

and /? = if and only if 

Tr(C*) = V * G iJ"(X, n\{D)). 

To make explicit the complex numbers Tr(C^), we use the Dolbeault 9-resolution 
of n\. We denote by P^'"^^ the sheaf of germs of (p, g)-currents on X. Let he a 
germ of meromorphic g-form aX, p £ X . We recall for convenience that the principal 
value current [-0] G 'D^x'p ^'^'^ residue current d[%lj\ G T^x'p'^ have the Cauchy 
integral representations 

([0],0) :=lim-^ / 

e^O 2nT Jun{\u\>e} 
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and 

e^O 2nT Jun{\u\=e} 

where 9 is some test-form with appropriate bidegree and u = is a local equation 
for the polar divisor of -0 in a small neighborhood U oi p (see [30j for instance). 
The Dolbeault resolution of is given by the exact complex of sheaves 

172, p(2,0) ^(2,1) p(2,2) ^ ^ 

This sequence breaks in the two short exact sequences 

Q^nl-.V^^^'^ ^Z'i''^ and 0-z(,^-^)^I?(,^'^)^I?(,^'^)^0, 

where -Zj^'^'' C I?^'^'' is the subsheaf of 9-closed (2, l)-currents. Since the sheaves 
I?^''^'* are fine, we obtain the two coboundary isomorphisms 

dH^(X,V)^' ') 

For any e H^{X,V?x{D)), there thus exists a global current T.^ G 2?^'^^) 
whose class [T.^] modulo 9 is unique solution to C* = '^2 ° '^ili^*]); and we have 
equality 

Tr(C*) = (T^,l). 
To make explicit Tip, we need the following lemma. 

Lemma 2. 1. There exists a Cartier divisor 7 on an open neighborhood B of \D\ 
which restricts to j on D. 

2. There exists a Cartier divisor E on X which restricts to T on Dred o,nd there 
is an isomorphism C ~ Ox (E) . 

Proof. 1. By mean of a partition of unity, we can construct a C°° function u on X 
which vanishes exactly on |£>|, giving an open tubular neighborhood of D 

B,^{\u\<e}. 

Let / = {fu} be a meromorphic section of C with Cartier divisor 7, and consider 
some local meromorphic liftings fij. We can choose V a sufficiently fine covering of 
B, so that fu e Ox{Uy for t/ n |r| = and so that n F n |r| = for distinct 
U,V ^ V. If now U intersects |r|, we can choose e' < e small enough so that fu 
has no poles nor zeroes on U O V O B^i. In such a way, fu/fv is invertible on 
U nV n B^r, giving a Cartier divisor 7 on B^i which restricts to 7 on D. 

2. Let !F — Ox{F) be some very ample line bundle on X with F an effective 
divisor which intersects properly |-Dred|- The isomorphisms 

combined with the structural sequence for Dred gives the exact sequence 

0^ C(^OxinF - Dred) -^C® Ox{nF) Od,.,A^ + nFi.) ^ 
for any integer n, where Tp — F ■ Dred- For n big enough, we have 

H\C®Ox{nF -Dred)) = ^ 

by Serre Vanishing Theorem so that F + nVp lifts to some Cartier divisor G on 
X . Then, the Cartier divisor E ~ G — nF restricts to F on Dred- By construction 
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Ox[E) and C ~ Od{i) have the same restriction to Dred = dX, and Proposition 
1 impHes £ ~ Ox(£^). □ 

Let _B,7 and E as in Lemma 2, with V the associated open covering of B. Since 
Ob{j) (8) Ob{~E) restricts to Cq on D, we can choose the hftings guv = "ijy/my, 
where m = {mu}ij^\> is the global meromorphic section of Ob{^) <E) Ob(— i?) with 
Cartier divisor 

div(m) = 7 - £^15 

on B. For _B small enough, rnmi/f^y G Og{U n y) for distinct U,V ^ V and there 
exists a logarithmic determination 



oo ^ 

\og{mu\unv) ■^y2~(^^ "^c/|;7ny)"- 



n=l 

Since div(m) ■ Dred = 0, the restriction of m to Dred is a constant, that we can 
suppose to be 1. There thus exists no £ N such that the meromorphic function 
(1 — m[/)" vanishes on the divisor D DU for all n > uq and all U E V. 

Let 4' e fl^{D)). We can multiply a principal value current (or a residue 

current) with a residue current as soon as their supports have a proper intersection. 
Then, by the Duality Theorem [3D], we have 

[{l-mu)"]d[^u] =0 Vn>no 
so that the principal value currents 



n 

n=l 

satisfy 

(6) (5c/ ~ Sv)d[^unv] = [log(5c/ny)]a[^'c/ny] 

for all U,V eV. We obtain in such a way a global (2, 2)-current Tb e r{B, P^^'^^) 
on i?, locally given by 

for all [/ G V. Since ^ is holomorphic outside \D\, we have Tjj = when [/nlDj = 
and we can extend Tg by zero to a global current T on X. Using (6) and the 
definition of the coboundary maps, we check easily that T is solution to 

62 o 6i{T) = C*. 

Let now V' be some rational 1-form on X so that d-ipixo — "^{Xo- For all [/ £ V, 
we define the local currents 



Ru := f dSu - \^] ) A BUu] and T(j = d\^] A BU 
\ I mu }/'-■' I mu J 



dmi 

mu 



From (6) we deduce equalities Ru\unv ~ Rv\unV: and the Ru's define a global 
(2, l)-current Rb on B. By assumption, mu/my is invertible on [/ n so that 
'^u\unv ~ -^viany Siting a global (2, 2)-current on B. Both currents i?B and 
are supported on iDj C B and can be extended by zero to global currents R 
and T' onX. 
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Since Tu has support a finite set, Stokes Theorem gives equalities 

{Tu,l) = {dSuAd[^u]A) 
= {dSu Add[^u]-l) 

Thus, (T, 1) = (T' + Br, 1) = (T', 1). Since div(m) = 7 - the Lelong-Poincare 
equation gives equahty 

(T', 1) = (T^j, 1) = ([7] A d[^lj]^B, 1) - {[E]iB A a[V^]|B, 1), 

where [■] designs here the integration current associated to analytic cycles. The cur- 
rent A9[V']|b is supported on the compact subset iFj C B, and the integration 
current is 9-closed. We deduce 

{[E]^B A 1) = {[E] A Biij], 1) = {dm A m, 1) - 0. 

The analytic set i? n 7 is a disjoint union of analytic cycles Jp — {fp = 0}. If V'p is 
the germ of at p in the chosen local coordinates, we obtain finally 

Tr(C*) - 5]([7p]Aa[V'p],l) 
pe|r| 

pe|r| ■'P pe\r\ 

and it follows that (2) is equivalent to that C — Od{i) extends to £ £ Pic(X). 
Note that the previous expression only depends on the hue bundle Od{i) and on 
^' by construction. If C extends to £, we can use Vanishing Serre Theorem as in 
Lemma 2 (with (I?, 7) instead of (DrediT)) and show that C ~ Ox{E) for some 
Cartier divisor E on X which restricts to 7 on D. The lifting bundle C being unique 
(up to isomorphism), the osculating divisor E is unique up to rational equivalence. 
This ends the proof of the first point. 

If conditions (2) hold, then C ~ Od{E) for some Cartier divisor E on X. By 
tensoring the structural sequence of D with C — Ox{E), we obtain the short exact 
sequence 

0-,Ox{E~D)^ Ox{E) ^ Od{E) ^ 0. 

If 7 is effective and H^{X, Ox{E - D)) = 0, the global section / e H°{D, Od{E)) 
with zero divisor 7 automatically lifts to / G H^{X,Ox{E)) and C :— divo(/) is 
an effective osculating divisor for (D,7). This ends the proof of Theorem 1. □ 

2.5. An explicit formula in the toric case. We show now that we can make 
explicit the conditions (2) in the case of a toric variety X. We refer the reader to 
[13] and [TT] for an introduction to toric geometry. 

2.5.1. Preliminaries. Suppose that X is a toric surface containing Xo as an union 
of orbits. Then, X has toric divisors Z^o, • • ■ , ^r+i where 

dX ^ Di + --- + Dr 

is the boundary of X and Dq and Dr+i are the Zariski closure of the 1-dimensional 
orbits of Xo ~ C^. 
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Let E be the fan of X. We denote by pi G E the ray associated to the toric divisor 
Di. Since E is regular, we can order the I?i's in such a way that the generators r]i 
of the monoids pi O 1? satisfy 

det (77^,77^+1) = 1 

(with convention r\r-\-i = ?7o)- We denote by Vi the afhne toric chart associated to 
the two-dimensional cone PiK"*" ® pi+iK"*". Thus, 

C/j = SpecC[a;j, ~ C^, 

where torus coordinates t — (ii,t2) and affine coordinates {xi,yi) are related by 
relations 

for all m = {7121,1712} E 1? , where — t^^t^^ (see [H]). With these conventions, 
the toric divisors Di and -Di+i have respective affine equations 

Di\Ui = i^i = 0} and A+i|;7i = {^i = 0} 
and \Dj | n C/^ = for i, « + 1. 

We can associate to any toric divisor E = X)i=o ^i^^i a polytope 
Pe = {m e {7n,7i,)+ei > 0, i = 0, . . . , r + 1}, 

where (■,■) designs the usual scalar product in K.^. Each Laurent monomial i"* 
defines a rational function on X with divisor 

r+l 
i=0 

and there is a natural isomorphism 

(7) H'iX,Ox{E)) ^ C • i'", 

where Ox (E) is the sheaf of meromorphic functions with polar locus bounded by 
E. Similary, the rational form defined by 

dtiAdt2 
^m|(c*)2-r — 

has divisor div(4'm) = div(t™) — Kx where —Kx := Dq + ■ ■ ■ + Dr+i is an anti- 
canonical divisor, and there is equality 

(8) H°iX,nj,iE))= C • 

mePE+Kx 

2.5.2. Explicit osculating criterions. Let (I?, 7) be an osculation data on dX. We 
note r = 7 ■ dX and Vi — •y ■ Di. By (8), the conditions (2) of Theorem 1 are 
equivalent to that 

(9) ^(7,^'„)p = VmePD+K^nZ^. 
pe|r| 

We give here an explicit formula for the residues (7, when the following hy- 

pothesis holds: 

(Hi) The zero-cycle F is reduced a7id does 7iot contain any torus fixed points. 
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In such a case, the germs 7p are smooth, irreducible, and have transversal intersec- 
tion with dX. For each p G |r|, there is an unique component Di containing p and 
we can choose the Weiertrass equation 

fp{xi,yi) ^Vi- <l>p{xi) 

for the germ 7p in the affine coordinates of Ui. The analytic function (jjp G C{xi} 
is well-defined modulo (a;^'"*"^) and does not vanish at 0. We obtain the following 

Proposition 1. Suppose that T satisfies (Hi) and let m G Pd+Kx ^ ^ ■ Then, 

1 a-('",'7.) 

<^'*"'-^ (-c..,.))!a.r'-4 CT)'"' 

for all p G \Ti\ and all i — 1, . . . ,r. We use here conventions (f>p/k :— log((/)p) for 
k^O and 8°- / dx" {(p^ / k) for a < 0. 

Since ^ Pd+Kxi equality (771,77^+1) = implies {m,rji) ^ and previous 
expression is well-defined since it depends only on the derivatives and logarithmic 
derivatives of 4>p evaluated at (recall that (/'p(O) ^ 0). Moreover, —{m,rii) < ki 
and (10) only depends on (f)p modulo (a;(^*^^) 

Example 2 (Wood's theorem). We keep the same hypothesis and notations of Ex- 
ample 1, with now an osculating order k > 2. Thus D = {k+ 1)L. Hypothesis {Hi) 
holds and by Proposition 1 , there exists an osculating divisor for {D, 7) if and only 
if 



(11) 9^^^(E i-)(0) = 0^ V7^G(N*)^ 7771 + 7772 <fc. 

p6|r' 



Since H^{Op2{d — — 1)) = it follows that (11) is also sufficient for the existence 
of an osculating curve. When k < d+1, we recover a theorem of Wood, 34J. 



2.5.3. Proof of Proposition 1. Let 771 G Pjj^Kx- I'^ particular (777,770) > 1 and 
we can suppose that 7712 ^ 0. The form is holomorphic in Xq and ^',n|Xo = 
d{4'm\Xo) where 4'm is the rational 1-form on X determined by its restriction 

dti 

^m|(c*)2 - —-r- 

7772 tl 

to the torus. Let p G iFij. We compute the associated residue in the chart Ui. If 
we let (ei, 62) be the canonical basis of R^, we obtain 

, {•m,rii) (m,j)i+i) 1 r/ i dXi . , dy 

Vm.\Ui=xl 'yl ^' — [(ei,7y,)— (ei,77,+i) — 



"72 Xi yi 



If {m,r]i) > 0, the form ijjm is holomorphic at p G iFij and reSp(y^ A -0™) 
Suppose now that {m,r]i) < 0. The Cauchy integral representation for Grothendieck 
residues of meromorphic forms does not depend on e for Ci small enough, thanks to 



= 0. 
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Stokes Theorem. Thus, for sufficiently smaU Ci's, we obtain 



reso 



x] "yl ■ -pA — 
Jp yi 



1 



(2i7r)2 
1 



\xi\=ei,\yi->t>p\=t2 



_.{'"^'?i>,,(™.')i+i> 



hp ^ dy. 



yi - 4'p yi 



(m,7)i-|-i) — 1 dyi ^ (mjii) 



yi - (pp 



1 



(-(m,77,)-l)! 



5 7 I (m,7?,+i) 



where the two last equalities are application of the Cauchy formula. In the same 
way, but simpler, we find 



reSr, 



(m,r?v) (m,77^ + i) dfp dXi 

Xi yl ^ A 



-1 



dx. 



'{ni^rji) \^P 



the minus sign coming from the chosen ordering of the numerator's factors. We 
deduce 



reso 



(f A*.)^ 



c 



1 



J Xi=0 



1112 {-(m^rii))! 
where 

C = (ei,7yi)(m,77i+i) - (ei, 77^+1 ) (to, 77i) = 7712 det(77i, 77^+1) = TO2. 
This ends the proof of Proposition 1. □ 

If two algebraic curves of fixed degree osculate each other on a finite subset 
with sufficiently big contact orders, they necessarily have a common component. 
We show in the next section that this basic fact permits to apply Theorem 1 and 
Proposition 1 to the Computer Algebra problem of polynomial factorization. 



3. Application to polynomial factorization. 

This section is devoted to develop a new algorithm to compute the absolute 
factorization of a bivariate polynomial /. The method we introduce is comparable to 
a toric version of the Hensel lifting. It is in the spirit of the Abu Salem-Gao-Lauder 
algorithm [T], but uses vanishing-sums as in the Galligo-Rupprecht algorithm [O] . 
Our main contribution is that the underlying algorithmic complexity depends now 
on the Newton polytope Nf instead of the degree deg{f). We prove our main 
results Theorem 2 and Proposition 2 in Subsections 3.2 and 3.3. We describe the 
algorithm in Subsection 3.4 and we comment and compare it with related results in 
Subsection 3.5. We discuss non-toric singularities in Subsection 3.6 and we conclude 
in the last Subsection 3.7. 
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3.1. Preliminaries and notations. Let / G K[ti,t2] be a bivariate polynomial 
with coefficients in some subfield K C C By absolute factorization of /, we mean 
its irreducible decomposition in the ring C[ti,t2], computed by floatting calculous 
with a given precision. 

Suppose that / has the monomial expansion 

We denote by Nf the Newton polytope of /, convex hull of the exponents to £ 
for which c„j 7^ 0. We recall that by a theorem of Ostrowski [28], we have relations 

for any polynomials /i, /2, where + designs here the Minkowski sum. 
We assume that the following hypothesis holds. 

{H2) The Newton polytope of f contains the origin. 

This equivalent to that /(0,0) ^ 0. 

Remark 1. We can always find s ^J? so that the Laurent polynomial ft^" becomes 
a polynomial f G C[ui, U2] with /'(O, 0) 7^ after some monomial change of affine 
coordinates (ti,t2) 1— > (wi,U2). This transformation preserves the sparse structure 
of f so that hypothesis {H2) is not restrictive for our purpose. 

An exterior facet of N f is a one-dimensional face F of Nf whose normal inward 
primitive vector has at least a negative coordinate. The associated normal ray of 
Nf is called an exterior ray. The associated facet polynomial fp of / is defined to 
be 

/f — ^ Cmt™". 

The facet polynomials have a one dimensional Newton polytope and become uni- 
variate polynomials after a monomial change of coordinates. 

We construct now a toric completion of associated to the exterior facets of 
Nf. To this aim, we consider the complete simplicial fan of whose rays are 

= {exterior rays of Nf} U {(0, 1)M+, (1, 0)M+}. 

The toric surface Xf = X^^ is a simplicial toric completion of — SpecC[ti,t2] 
whose boundary dXf = Xf\C'^ is the sum of the irreducible toric divisors associated 
to the exterior rays of Nf. The intersection of the Zariski closure C/ C Xf of the 
affine curve {/ = 0} C with the irreducible components of the boundary dXf 
corresponds to the non trivial roots of the exterior facet polynomial of / (see [23] 
for instance). In particular, the curve Cf does not contain any torus fixed points 
(it's clear for those contained in dXf and the remaining torus fixed point of Xq 
does not belong to Cf by {H2)). 

We say that / satisfies hypothesis (Hi) of Subsection 2.5 when the zero-cycle 
Tf = Cf ■ dXf does. This corresponds to the case of exterior facet polynomials of 
/ with a square free absolute decomposition in C[i,t^^]. 
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The boundary of Xf contains the singular locus S := Sing{Xf) of Xf and is 
generally not a normal crossing divisor. To this aim, we consider a toric desingu- 
larization 

X^Xf 

of Xf, whose exceptional divisor E satisfies 5 = 7r(|i?|). The smooth surface X is 
now a projective toric C^-completion with a toric normal crossing boundary 

dX = Di + --- + Dr 

whose support contains \E\. Since Cfr\S = 0, the total transform C = 7r"i(C/) 
of Cf has a proper intersection with dX and the irreducible decomposition of C 
corresponds to the absolute factorization of /. Note that C is not necessarily 
reduced. 

As in Subsection 2.5, we denote by po, . . . , Pr+i the rays of the fan S of X, and 
by r/o, • • • , Vr+i their primitive vectors. So r]o = (0, 1) and r]r+i = (1, 0). The curve 
C has a positive intersection C ■ Di > when pi e S(l) is an exterior ray of Nf, 
and C • A = when pi G E(l) \ (that is when 7r(|A|) C S). 

The following lemma will be usefuU: 

Lemma 3. There is an unique effective divisor D linearly equivalent to C + dX 
with support \dX\ and 

Pd+kx nz^ =iV/n(N*)2. 

Proof. The polynomial / extends to a rational function on X whose polar divisor 
diV(x)(/) is supported by \dX\. Thus, the divisor 

£»:=diVoo(/) + aX 

is effective, with support \dX\, and rationally equivalent to C + dX. If D' is an 
other candidate, then D — D' = div(/i) for some h G C{X) with no poles nor zeroes 
outside \dX\. Since X \ \dX\ ~ C^, the rational function h is necessarily constant 
and D = D'. By using afBne charts, we check easily that the polar divisor of / is 
equal to 

r 

divoo(/) = V fciA, ki = - min (m,r?j). 

'—^ m£Nf 
i=l ' 

By (7), it follows that Nf is contained in the polytope Pg of G = divoo(/) and 
intersects each face of Pq ■ Since the fan of X refines noth normal fans of G and Nf 
(we use here hypothesis {H2)), it follows that Pq = Nf. Then, equality D + Kx = 
G — Do — Dr+i implies that 

Pd+k^ = {m e Pg, (m,r?o) > 1, (m,r/^+i) > 1} = Pg n (N*)^. 

□ 



3.2. Computing the Newton polytopes of the absolute factors. We show 
here how to recover the Newton polytopes of the irreducible absolute factors of Nf. 

We recall that the Newton polytope of a factor of / is necessarily a Minkowski 
summand of Nf. For a general polytope Q, we denote by Q''*'' the set oi m G Q for 
which the scalar product {m,r]i) is minimal. 

We obtain the following 
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Theorem 2. 1. Let f G K[ti,t2] which satisfies (i?2)- Let tt, X and C be as 
before and denote by 7 the restriction of C to the divisor D of Lemma 3. Let Q be 
a lattice Minkowski summand of Nf . 

There exists an absolute factor q of f with Newton polytope Q if and only if there 
exists j' < J an effective Cartier divisor on D with 

(12) deg(7' • A) = Card(g(*) n Z^) - 1, i = l,...,r 

and such that (2) holds for the osculation data [D,^') on the boundary of X . The 
polynomial q can be recovered from 7' . 

Note that Q(*) is reduced to a point and Card{Q''^^ fl Z^) — 1 = as soon as pi 
is not an exterior ray of Nf. 

Proof. A Minkowski sum decomposition Nf = P + Q corresponds to a line bundle 
decompostion 

where G, Gp and Gq are the polar divisors of the rational functions determined by 
polynomials with respective polytope Nf, P and Q. Note that G = Gp + Gq. Since 
S refines the normal fan of each polytope, we have equalities Nf = Pg, Q — Pgq 
and P = Pqp (see Lemma 3) and all involved line bundles are globally generated 
over X (see [I4 ). 

A factor of / with Newton polytope Q corresponds to a divisor G' < G lying 
in the complete linear system \Ox{Gq)\. It defines a Cartier divisor 7' = Cj'^ < 
G\D = 7 on _D. Thus, conditions (2) hold for the osculation data {D,j') and 

deg(7' • A) = deg{G' ■ A) = deg{Ox{GQ)\D,) 

for alH = 1, . . . ,r. Since Ox{Gq) is globally generated over X, toric intersection 
theory gives equalities 

deg(Ox (Gq) I i5 J - Card(g« n Z^) - 1 

for all i = 1, . . . , r. This shows necessity of conditions. 

Suppose now that {D, 7') satisfies (12) and (2) and let F be an osculating divisor 
asociated to that data. We deduce equalities 

deg(Ox(-F)|A) = deg(7' ■ A) = deg{OxiGQ)\Dj 

for all i = 1, . . . , r. Thus, there is an isomorphism 

Oax{F) ~ OaxiGQ) 

and OxiF) ~ Ox(Gq) by Corollary 1 of Subsection 2.4. Finally, (2) and (12) hold 
for (£>, 7') if and only if Od{j') — Od{Gq). Let us show that there exists in such a 
case an effective osculating divisor. The problem is that H^{X, Ox{Gq — D)) does 
not necessarily vanish. To avoid this difficulty, we introduce the effective divisor 

Dq := Gq + dX < G + dX = D. 
Then OD^iGg) Oi5j7pJ and H\X,Ox{Gq - Dq)) ^ H\X,Ox{-dX)). 
Since dX is reduced and connected, the restriction map H°{Ox) H"{Odx) is an 
isomorphism and we deduce equality H^{X, Ox{—dX)) — 0. By Theorem 1, there 
thus exists G' e\Ox{GQ)\ with 

^/ / 
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Such a curve has proper intersection with the remaining toric divisors Dq and Dr+i 
(otherwise the support of 7' would contain a torus fixed point). We deduce that 
there exists a unique Cartier divisor 70 on D + Dq + Dr+i so that 

lo\D = 7' and 7o|Do+D,+i = C ■ Dq + C ■ Dr+i 
We have equality D + Dq + D^+i = G — Kx and we obtain 

H\X,Ox{GQ-D-Da-Dr+i) ~ H\X,Ox{Gq^G + Kx) 

~ H\X,Ox{G-GQ)j 
~ H\X,Ox{Gp))^0. 

The second isomorphism is Serre Duality, and the last equality comes from the fact 
that the of a globally generated line bundle on a toric surface vanishes (see 
[14j). Last equality combined with the short exact sequence 

O^Ox{Gq-D-Do- Dr+i) ^ Ox{Gq) ^ Od+d„+d^^Aio) 0, 

implies that there exists G" G \Ox{Gq)\ which restricts to 70 on Z? + Dq + Dr+i- 
A fortiori, G" restricts to 7' on D. Since Gq — Dq = —dX < 0, we have 
H°{X,Ox{Gq - Dq)) = 0. Thus G' ^ G" , the two curves having the same 
restriction to Dq. 

There remains to show that C" < G. Suppose that there exists an irreducible 
component Gq of C" with proper intersection with G. Then, we can compute the 
intersection multiplicity of Gq and G at any p E X. It is given by 

mult„(Co, G) = dime 

(/p:5p) 

where fp and gp are respective local equations for G and Go at p. By assumption, 

Go\D < G\D 

so that the class of gp in Od,p divides that of fp. We deduce that 

fp — Upgp -\- Vphp 

for some holomorphic functions Up, Vp, where hp is a local equation for D. Thus, 
we obtain inequality 

dime . > dime , '^^'^ , = multp(Go, D), 

Up,9p) {gp,hp) 

for all p in the support of D. Finally, 

deg(Ox(Go))|c > multp(Go,C)> ^ multp{Go, D) = deg{Ox{Go))iD- 

p£\D\ pe\D\ 

Since D is rationally equivalent to G + dX, this implies that Gq has a negative 
intersection with the boundary. This leads to a contradiction. Thus Go is neces- 
sarily an irreducible component of G. Finally we have G' < G, corresponding to 
an absolute factor q oi f with Newton polytope Q. 

Since Gq - Dq < 0, the restriction H°{Ox{Gq)) H°{Odq{Gq)) is injective. 
Thus G' e \Ox{Gq)\ can be uniquely recovered from its restriction 7p^^ to Dq 
(and a fortiori from 7'). □ 
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Remark 2. // we only assume that Nf intersects both coordinate axes, the oscu- 
lating criterions still permit to detect an absolute factor q of f associated to the 
choice of a summand Q of Nf . Nevertheless, although the Newton polytope Nq of 
q still has the same exterior facets of Q, we can not conclude that Nq = Q. 

3.3. Computing the absolute factors of /. Suppose that 7' < 7 satisfies con- 
ditions of Theorem 2 for a Minkowski summand Q of Nf. We give here an efficient 
way to compute the corresponding absolute factor g of /. The polynomial q admits 
the monomial C-expansion 

q{t)= Yl 

meQnN2 

and we look for the homogeneous class [a] G pCard(QnN of a = (am)meQnN2- 

We denote by T' = 7' • dX and by T^ = T' • A- Let us fix p G |r^|. If / satisfies 
the hypothesis {Hi ) , the germ of C at p has the Weirstrass equation 

Ui - (kpixi) = 0, 

and we can define 

ap{u,v) :— — — ^(0) if u > 0, apiu,v) := if u < 0, 
ax" 

for all integers u, v (recall that 0p(O) ^ 0). For fc G N and m G Z^, we then define 
the complex number 

/3p{k,m) := ap{k - {m,r]i) - a, {m,r]i+i) + Ci+i), 

where := — minmgQ(rn, 77^). We denote by Yol{Q) the euclidean volume of Q in 
M^. We obtain the following 

Proposition 2. Suppose that f satisfies (Hi) and {H2). The homogeneous vector 
[a] is uniquely determined by conditions 

for all p G |r^|, all < k < e.; and all i — 1, . . . ,r. The underlying linear system 
{S^>) contains 2 Vol((3) + deg(r') equations of Card((5 H Z^) unknowns and only 
depends on the restriction of j' to the divisor Gq + dX . 



Proof. By (11] , the curve C of q has affine polynomial equation qi — Q 'm the chart 
Ui = Spec C[xi, T/i], where 

/ \ \ " (m,rii)+ei (m,77i + i)+ei+i 

qt{x,,y^)= a.^x\ y\ +. 

rriGQnZ^ 

Let p G |r^|. Since the germ of C" at p is contained in that of C it follows that 

qi{x,i,(i)p{xi)) = 

It's easy to show that the fc*''-derivative of this expression evaluated at is equal 
to c X X]meQnz2 ctmPpik^m) for a non zero scalar c. Thus the coefficients of [a] 
determine a non trivial solution of {S^'). 
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If a = (am)meQnz2 is an other non trivial solution of {Sy), the polynomial 
Qit) = Emegnz^ Qm^" satisfies 



d' 



k 



= 



for all p e |r^|, all < fc < Bi and alH = 1, . . . , r. By the duality theorem, this 
forces the germ of §i at p to belong to the ideal {yi — (j)p{xi), x^'^^). It follows that 
the curve C C X oi q satisfies 

Since Ng C Q — Ng, it follows that deg{C ■ Di) < deg{C' ■ Di) for all i. Combined 
with previous inequality, this forces equality 

We have chosen the e^'s in order to that 

r 

Y,{e, + l)D, = GQ + dX, 

i=l 

and both curves C and C belong to the linear system \Ox{Gq)\. The restric- 
tion map H^[X,Ox{Gq)) — > H^{X,OGQ+dx{GQ)) being injective, the previous 
equality forces C — C. It follows that q = cq for some c S C* and [a] = [a]. 

By construction, the linear system {S^') only depends on the restriction of 7' to 
the divisor Gq + dX < G + dX = D. It contains precisely 

r r 

^ ^ (e, + 1) = ^(e, + 1) deg(C" • A) - deg(C' • Gq) + deg(r') 

i=lpG|r'J i=l 

equations and deg(C" • Gq) = deg(C" ■ C") = 2 Vol((3) (see [IJ for instance). □ 



Remark 3. Each linear equation in Proposition 2 involves a reduced number of 
unknowns and the linear system {Sj') has a very particular sparse structure. For 
instance, letting k — 0, we obtain for each i ~ 1, . . . , r the linear subsystem 

It has Card{Q^'^^ H 7?) — 1 equations with Card{Q^^^ H 7?) unknowns and permits 
to determine the coefficients of the i^^ exterior facet polynomial of q. For k — 1, 
we deduce relations on the coefficients {am, (^,77^) + Ci = 1}. For a general k, we 
deduce relations on the coefficients {am, (m,rii) + = k}. In general, there are 
no non trivial solutions to (S-yi). The vector subspace of solutions has dimension 1 
precisely when 7' obeys to the osculation conditions (2) . 



OSCULATION. 



21 



3.4. A sparse vanishing-sums algorithm. We describe here a sparse vanishing- 
sums algorithm associated to Theorem 2 and Proposition 2. When we say compute, 
we mean compute by using floatting calculous with a given precision. When we say 
test a vanishing-sum, we mean test an < e-sum for an arbitrary small e > 0. 

Input: A polynomial / G K[ti,t2] which satisfies /(0,0) / and with square free 
exterior facet polynomials. 

Output : The irreducible factorization of / over C. 

Step 1. Compute the Minkowski summands of Nf. Use for instance the algorithm 
presented in [TO. If Nf is irreducible, so is /. Otherwise go to step 2. 

Step 2. Compute the fan of X . The fan E is obtained from "Ef hy adding some 
rays in the singular two-dimensional cones of . Such a fan can be obtained in a 
canonical way by using an Euclidean algorithm, or by computing some Hirzebruch- 
Jung continued fraction (see [14 ). 

Step 3. Compute the osculating data {D,"f). By Lemma 3, we have equality 

r 

D = y^(ki + fcj := - min (m, 77.,). 

^ — ' m^Nf 
i=l ' 

The Cartier divisor 7 on D is obtained by computing the family of implicit functions 
{(j)p, p G iTil} to order fc^, for i = 1, . . . ,r. Note that Fj = when pi is not an 
exterior ray of Nf. 

Step 4- Compute the Newton polytopes of the absolute irreducible factors. Proposi- 
tion 1 combined with Theorem 2 gives an efhcient way to compute the decomposi- 
tions 

7 = 7i-H.---h7,, Nf = Qi + --- + Q, 
of 7 and Nf associated to the irreducible absolute decomposition f — qi ■ ■ ■ Qs of f. 

Step 5. Compute the irreducible absolute factors of f. Use Proposition 2. The 
coefficients of the linear systems (S'-yJ, i — 1, . . . , s are obtained from the residues 
(7, *m)p, m G QiC] (N*)^ already computed in step 4. 

The numerical part of the algorithm reduces to the computation of the roots 
of the univariate exterior facet polynomials of /. Then it detects the absolute 
factorization of / with a probability which increases after each positive vanishing- 
tests, up to obtain the adequat decomposition of Nf. Finally, we compute the 
factors by solving some linear systems. Roughly speaking, we obtain here a toric 
version of the Hensel lifting (see Remark 3). A comparable algorithm has been 
obtained by Abu Salem-Gao-Lauderin in [1], by using combinatorial tools. 

If we can decide formally if a sum vanishes, there is no chance of failure in step 
4 and the algorithm is deterministic. If we test the osculating criterions (9) for a 
generic linear combination of the involved m's, the Newton polytope decomposition 
is valid only with probability one, in the vain of the Galligo-Rupprecht and Elkadi- 
Galligo-Weimann algorithms or |12j . 

As in [13] and [12], our algorithm necessarily uses < e-sum tests imposed by 
floatting calculous and numerical approximation. Thus, it can happen that the 
Newton decomposition of Nf in step 4 does not correspond to the absolute decom- 
position of / (and that for any choice of e > 0) and there is a chance of failure 
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of the algorithm. Nevertheless, in the important case a polynomial / defined and 
irreducible over Q the authors in [7J show that we can recover the exact factor- 
ization with formal coefficients in a finite extension of K from a sufficiently fine 
approximate factorization. This problem will be explored in a further work. 

3.5. Comparison with related results. We compare here our algorithm with 
that of Galligo-Rupprecht [13] and with that of Elkadi-Galligo-Weimann [12] . 

3.5.1. Comparison with the Galligo-Rupprecht GR-algorithm. In |13j . the authors 
perform a generic change of affine coordinates and then compute the factorization 
of/ 

d 

1=1 

in C{ti}[<2] modulo {t\). They detect the factors of / with probability one by 
testing the Reiss relation (3) on subfamilies T C {(f)i, i — 1, . . . ,d}. Then they lift 
and compute the candidate factors by using Hensel lemma. Let us compare with 
our approach. 

The recombination number. We define the recombination number M{f) to be the 
maximal number of choices 7' < 7 necessary to detect the absolute factorization of 
/ in step 4. By (12), it depends on the geometry of Nf and is subject to constraints 
given by the possible Minkowski-sum decompositions of Nf. For instance, if / is 
irreducible over its irreducible absolute factors necessarily have the same polytope 
(see [8] for instance) and N{f) decreases drastically. 

We denote by M{f) the recombination number of the GR-algorithm, that is the 
number of possible choices for the families J- when taking in account restrictions 
imposed by the possible Minkowski-sums decompositions of Nf. 

We have the following improvement: 

Lemma 4. Suppose that f is irreducible over K and satisfies hypothesis (Hi) and 
{H2). Letd^degif). Then 

M{f) < M{f) < 2'^ 

with equivalence 

M{f) = M{f) ^ iV/ = Conv{(0,0),(d,0),(0,d)}, 
where Conv designs the convex hull. 

Proof, let n be the biggest integer so that Nf — nQo for a lattice polytope Qo- 
Thus d — nl where I is the total degree of a polynomial with polytope Qo- By [8], 
an irreducible absolute factor of / has Newton polytope fcQo for some integer k 
which divides n and the GR-algorithm looks for an irreducible factor of degree kl. 
This gives the recombination number 

k\n 

where Cj is the usual number of combinations. We clearly have M{f) < 2'^. 
Equality M{f) = 2'^ holds if and only if / is a dense polynomial which is not 
assumed to be irreducible over K. 
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Denote now by li, . . . ,lt the lattice lenght of the exterior facets of Qq. Then, 
restrictions (12) imposed to the possible choice of 7' < 7 induce equality 



k\n i—1 



We have both inequalities 



n^nh ^ ^n(/i+---+it) 1 ^n(/i + ---+/t) ^ (-ml 

^kU ^ ^k{l^+-+U) ^k{h + -+lt) - ^kl 

i=l 

for any k < n. First inequality is an equality if and only if there is only t ~ 1 
exterior facet, and the difference strictly increases with t. Second inequality follows 
from the fact that the sum 

If — 7111 + ■ ■ ■ + nit 

of the lattice lenght of the exterior facets of Nf is smaller or equal to d = deg(/). 
Moreover, we can convince that there is equality If = d if and only if the normal 
fan S/ of Nf is regular. This is of course exceptional, and If « 0? in general. 
Finally, we have M{f) < M{f) and equality holds if and only if Nf is regular 
with one exterior facet. Since G Nf by assumption, this is equivalent to that 
Ar/ = Conv{(0,0),(d,0),(0,d)}. □ 

Note that there are fast factorization algorithms over a number field K (see [4], 
[26]). The following example illustrates Lemma 4. 

Example 3. Suppose that / is irreducible over K and that Nf — uQq, where Qq is 
the "undivisible" lattice polytope 

Qq = Conv{(0, 0), (a - 1, 0), (0, a -I), (a, a)}, 

with a > 2 and n prime. The lattice lenghts of the t = 2 exterior facets of Nf are 
both equal to n. Since n is prime, the underlying recombination number is equal 
to 

Mif) = 

In particular, is does not depend on a. The GR-algorithm considers / as a dense 
polynomial of degree deg{f) = 2na and looks for an irreducible absolute factor of 
degree 2a. The induced recombination number A4{f) is thus equal to 

and growths exponentially with a. 

The asymptotic estimation of Af{f) for a "generic polytope"0 Nf is a difficult 
problem. It is related to the estimation of the number of exterior facets, and to 
the estimation of "how singular" is the normal fan. Geometrically, these numbers 
correspond to the Picard numbers of Xf and of a minimal resolution X . 

Note that by Lemma 3 and Theorem 2, the maximal number of vanishing-sums 
to test is equal to 

Card{Nf n {Wf) x M{f). 



^The genericity has to be defined relatively to some invariant, as the cardinality of interior 
lattice points, or the volume for instance. 
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If we only want to detect absolute factors with probability 1, it's enough to test (9) 
on a generic linear combination of the involved m e Nf Hi?. In such a case, it's 
enough to test A/'(/) vanishing-sums. 

The numerical part. The numerical part of our algorithm reduces to the computa- 
tion of the roots of the exterior facet univariate polynomials. It's faster to factorize 
t univariate polynomials of degree li, . . . ,lt than an univariate polynomial of total 
degree d > li + ■ ■ ■ + It- Thus, the computation of C ■ dX is faster than with a 
generic line as soon as 7V(/) < M{f). 

The lifting step. Although we need to compute a reduced number implicit func- 
tions, we need in general a bigger precision on the 0p's than the upper bound 
deg(/) precision required with the classical Hensel lifting. In Example 3, we need 
for instance to compute the (pp^s up to order na^ >> 2na — deg(/). Morally, the 
more the recombination number is reduced, the more the required precision on the 
^p's increases. Thus, we should be very carefuU when comparing algorithmic com- 
plexity. Nevertheless, we gain on both sides in the important case of a polynomial 
of bidegree (a, 6). Formulas (9) and (10) shows that we need to compute a and b 
implicit functions up to respective maximal orders b and a, while the GR-algorithm 
computes a + b implicit functions up to a maximal order a + b. 

3.5.2. Comparison with the Elkadi-Galligo-Weimann EGW- algorithm. In [12j . the 
authors develop the sketch of an algorithm with the same recombination number 
that here by using the interpolation criterions obtained in |31| . Their approach 
necessites to compute numerically the intersection of C C X with a generic curve 
L in a very ample linear system which is "close enough" to \dX\. In other words, 
they pick a Newton polytope P whose normal fan is that of X and then solve 
a polynomial system / = ep + 1 = 0, where p has polytope Np — P and e is 
a small positive real number. When e goes to 0, the roots of the system go to 
the boundary of X and they deduce an asymptotic distribution of the zero-cycle 
C ■ L which traduces the polytopal information. As in the GR-algorithm, the 
generic choice of p in the EGW-algorithm permits to compute deg(C • L) implicit 
functions only up to order 2 in order to detect the absolute decomposition of / 
with probability one. Roughly speaking. Theorem 2 corresponds to the limit case 
e = 0. A great advantage is that we avoid the delicate problem of the "small 
enough e" choice and of the asymptotic distribution lecture. We compute only 
deg(C • dX) « deg(C ■ L) implicit functions but with a precision >> 2 and we 
detect the absolute decomposition of / deterministically. 

3.6. Using non toric information. Suppose that C is reduced, but with singu- 
larities along the boundary of X (so the exterior facet polynomials of / have square 
absolute factors). Thus F ~ C -dX is non reduced and the computations of residues 
is much more delicate. In particular, we can not use formula (10). 

Nevertheless, there exists in such a case a (non toric) smooth completion X of 
obtained from X by a serie of blow-ups, and so that the proper transform C C X of 
C has a transversal intersection with the boundary of X. By choosing an effective 
divisor D supported on dX and with sufficiently big multiplicities, we can then use 
Theorem 1 efficiently to decompose C and to recover the absolute factorization of 
/. The added exceptional divisors give new restrictions on the possible choices of 
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7' < 7 (where 7 = C|g) and the presence of singularities of C along dX finally 
turns out to be an opportunity to reduce the recombination number. There remains 
to find the best choice for D. This will be explored in a further work. 

3.7. Conclusion. Wc propose a new algorithm which computes the absolute fac- 
torization of a bivariatc polynomial by taking in account the geometry of the New- 
ton polytopc. For a sparse polynomial, this permits to reduce the recombination 
number when compared to the usual vanishing-sums algorithms. There remains to 
implement such an algorithm and to compute formally its complexity. What we'll 
remind here is the general idea that 

The more a curve is singular, the more it is easy to decompose. 
For a sparse polynomial /, a naive embedding of the curve of / in produces many 
toric singularities on the line at infinity, and we have shown here that osculation 
criterions in an appropiate toric resolution X permit to use this sparse information. 
Moreover, Theorem 1 permits in theory to profit also of the non toric singularities 
of C on the boundary of X. 

References 

[1] F. Abu Salem, S. Gao, A.G.B. Lauder, Factoring polynomials via polytopes, proc. of ISSAC 
(2004), pp. 411. 

[2] M. Andcrsson, Residue currents and ideal of holomorphic functions, Bull. Sci. math. (2004), 

pp. 481-512. 

[3] M. Avcndano, T. Krick, M. Sombra, Factoring bivariate sparse (lacunary) polynomials, J. 

Complexity 23 (2007), pp. 193-216. 
[4] K. Belabas, M. van Hoeij, J. Kluners, A. Steel, Factoring polynomials over global fields, 

Manuscrit (2004). 

[5] W.P. Barth, K. Hulek, C.A.M. Peters, A. Van De Ven, Compact Complex Surfaces, Springer- 
Vcrlag, Second Edition (2004). 

[6] A. Bostan, G. Lecerf, B. Salvy, E. Schost, B. Wiebelt, Complexity issues in bivariate poly- 
nomial factorization, proc. of ISSAC (2004), pp. 42-49. 

[7] G. Chcze, A. Galligo, From an approximate to an exact factorization, J. Symbolic Compu- 
tation, 41, no. 6 (2006), pp. 682-696. 

[8] G. Chczc, Absolute polynomial factorization in two variables and the knapsack proble, proc. 
of ISSAC (2004), pp. 87-94. 

[9] G. Cheze and G. Lecerf, Lifting and recombination techniques for absolute factorization, J. 
of Complexity 23, no. 3 (2007), pp. 380-420. 
[10] D. Cox, The homogeneous coordinate ring of a toric variety, J. Algebraic Geom. 4 , no. 1 
(1995), pp. 17-50. 

[11] V. Danilov, The geometry of toric varieties, Russian Math. Surveys 33 (1978), pp. 97-154. 
[12] M. Elkadi, A. Galligo, M. Wcimann, Towards Toric Absolute Factorization, ,1. Symb. Comp. 

(2009), doi:12.1216/j.jsc.2008.03.007. 
[13] A. Galligo, D. Rupprccht, Irreducible decomposition of curves, J. Symb. Comp., 33 (2002), 

pp. 661-677. 

[14] W. Fulton, Introduction to Toric Varieties, Annals of Math. Studies, Princeton University 
Press (1993). 

[15] S. Gao, Factoring multivariate polynomials via partial differential equation. Math. Comp. 

72, no 242 (2003), pp. 801-822. 
[16] S. Gao, A.G.B. Lauder, Decomposition of polytopes and polynomials. Disc, and Comp. Geom. 

6, nol (2001), pp. 89-124. 
[17] ,1. von zur Gatlicn, J. Gerhard, Modem computer algebra, Cambridge University Press, 1st 

edition (1999). 

[18] M.L. Green, Secant functions, the Reiss relation and its converse. Trans, of Amer. Math. 

Soc. 280, no.2 (1983), pp. 499-507. 
[19] P.A. Griffiths, Variations on a theorem of Abel, Inventiones Math. 35 (1976), pp. 321-390. 



26 



MARTIN WEIMANN 



[20] P.A. Griffiths, J. Harris, Principles of Algebraic Geometry, Pure and applied mathematics, 

Wiley-Interseiences (1978). 
[21] P.A. Griffiths, J. Harris, Residues and zero-cycles on algebraic varieties, Annals of Math. 

128 (1978), pp. 461-505. 
[22] R. Hartshorne, Algebraic Geometry, Springer- Verlag (1977). 

[23] G. Henkin, M. Passare, Abelian differentials on singular varieties and variation on a theorem 

of Lie-Griffiths, Inventiones math. 135, 297-328, 1999. 
[24] A.G. Khovansky, Newton polyhedra and toric varieties, Punct. Anal. Appl. 11 (1977), pp. 

56-67. 

[25] G. Leeerf, Sharp precision in Hensel lifting for bivariate polynomial factorization. Math. 

Comp. 75 (2006), pp. 921-933. 
[26] G. Leeerf, Improved dense multivariate polynomial factorization algorithms, J. Symb. Comp. 

(2007), to appear. 

[27] A.K. Lenstra, H.W. Lenstra, Jr. and L. Lovz, Factoring Polynomials with Rational Coeffi- 
cients, Math. Ann. 261, no.2 (1982), pp. 515-534. 

[28] A.M. Ostrowski, On multiplication and factorization of polynomials. Lexicographic orderings 
and extreme aggregates of terms, Aequationes Math. 13 (1975), pp. 201-228. 

[29] D. Rupprecht, Semi-numerical absolute factorization of polynomials with integer coefficients, 
J. SymboUc Computation 37 (2004), pp.557-574. 

[30] A. Tsikh, Multidimensional residues and their applications, Transl. Amer. Math. Soc 103, 
(1992). 

[31] M. Weimaim, Trace et calcul residuel: nouvelle version du theoreme d'Abel-inverse et formes 

abeliennes, Ann. de Toulouse, 16, no.2 (2007), pp. 397-424. 
[32] M. Weimann, An interpolation theorem in toric varieties, Ann. Inst. Fourier 58, no.4 (2008), 

pp. 1371-1381. 

[33] J. A. Wood, A simple criterion for an analytic hypersurface to be algebraic, Duke Math. J. 
51 (1984), pp. 235-237. 

[34] J. A. Wood, Osculation by Algebraic Hypersurfaces, J. Differential Geometry 18 (1983), pp. 
563-573. 

Departament Algebra i Geometria, Facultat de Matematiques, Universitat Barcelona 
Gran via 585, 08007 Barcelona. 

E-mail address: weimann23Sgmail.com 



